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^ ■ Abstract 

, We present a quenched lattice QCD calculation of the first few moments of the polarized and 

o , 

lO ' un-polarized structure functions of the nucleon. Our calculations are done using domain wall 

|: „,.„.,... 



X 



V ~ (2.4 fm)'^, and light quark masses down to about 1/4 the strange quark mass {m-n- ~ 400 
MeV). Values of the individual moments are found to be significantly larger than experiment, as 
in past lattice calculations, but interestingly the chiral symmetry of domain wall fermions allows 
for a precise determination of the ratio of the flavor non-singlet momentum fraction to the helicity 
distribution, {x)u~d/{x)Au-Aa^ which is in very good agreement with experiment. We discuss the 
implications of this result. Next, we show that the chiral symmetry of domain wall fermions is 
useful in eliminating mixing of power divergent lower dimensional operators with twist-3 operators. 



Finally, we compute the isovector tensor charge at renormalization scale // = 2 GeV in the MS 
scheme, (l)^^-^^ = 1.192(30), where the error is the statistical error only. 

PACS numbers: ll.15.Ha, 11.30.Rd, 12.38.Aw, 12.38.-t 12.38. Gc 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD) is the theory describing the strong interactions, and 
hence it is responsible for the properties of hadronic matter. Unhke Quantum Electrody- 
namics (QED), its non-perturbative nature, or strong coupling constant, makes it difficult to 
understand the low energy content of the theory. The lattice formulation of QCD provides 
both a non-perturbative way of defining the theory and a very powerful tool to calculate its 
properties. 

Deep inelastic scattering of leptons on nucleons has been the basic experimental tool in 
probing QCD [ll O, Is, 0, O, 0, Q, 0, 0, [l^ 11, These experiments have given rise to 
the parton model and, through extensive fits, have indirectly allowed the measurement of 
the parton distribution functions, the basic structural blueprint for hadrons. Connecting 
these experiments to the underlying theory of QCD is an important theoretical endeavor. 
During the last few years lattice computations have provided many interesting results for 



nucleon matrix elements 



i lattice computations nave provided many inl 

Q, 14, H 0, Q, E, 22, H Q 



in both quenched 



and full QCD. These calculations provide first-principles values for the moments of structure 
functions at leading twist. One of the major unresolved issues in these previous calculations 
is the approach to the chiral limit; computational limitations have restricted calculations 
to relatively large quark masses, introducing ambiguities in the extrapolation to the chiral 



limit. Furthermore, values calculated using the lattice regu 



arization have significantly over- 
25l | , leading to suggestions in the 



literature that strong suppression in the chiral limit is required to resolve the problem 



We address this question with a calculation using domain wall fermions 
Preliminary results have been given in 



17 



19 



28 



26|. 
291. 



The use of domain wall fermions allows 



us to examine the source of several systematic errors. Chiral symmetry at non-zero lattice 
spacing minimizes discretization errors, O(a^) in this case. Thus we work with relatively 
coarse lattice spacing, a ~ 0.15 fm, and therefore larger physical volume L ^ 2.4 fm. This 
means calculations with light quark masses will not suffer unduly large finite size corrections. 
In this study the lightest quark mass is roughly 1/4 of the strange quark mass, as light as has 
been used in nucleon structure calculations. Chiral symmetry makes the renormalization 
properties of operators simpler since there is less mixing with unwanted operators. Here 
operators are non-perturbatively renormalized, reducing a significant source of systematic 
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error. We have also chosen to use the DBW2 gauge as it substantially reduces the already 
small explicit chiral symmetry breaking for domain wall fermions with finite extra fifth 



The remainder of this paper is organized as follows. In Section O we briefiy recall the 
polarized and unpolarized structure functions of the nucleon arising from deep inelastic 
scattering and the operators that arise from their operator product expansions. The lattice 
transcription of operators and correlation functions are described in Section IIIII Perturba- 
tive and non-perturbative aspects of operator renormalization are discussed in Section IIVI 
Details of the numerical simulation are given in Section |3 Section Wl\ containing the pre- 
sentation and discussion of results, is the main part of the paper. We summarize the present 
study and comment on future calculations in Section rvill 

II. NUCLEON STRUCTURE FUNCTIONS 

The cross-section for deep inelastic scattering of leptons on a nucleon target is given by 
the square of the matrix element for an initial state lepton-proton pair to scatter to a final 
state of a lepton and hadrons. After summing over all possible final states, the square of the 
matrix element is computed using the optical theorem which relates the summed, squared, 
matrix element to the forward matrix element between nucleon states of the product of two 
electromagnetic currents. 



where L^*^ and W^^, are leptonic and hadronic tensors respectively, and q is the space-like 
four-momentum transfered to the nucleon by scattering off the electron through a virtual 
photon. The leptonic part is handled in perturbation theory since the QED coupling constant 
is small; the hadronic part, however must be treated non-perturbatively which is the focus 
of this paper. 

The hadronic tensor is conveniently split, PV^j, = W^^'^^ + W^'^'^\ The symmetric piece 
defines the unpolarized, or spin-average structure functions Fi and F2 (and F3 if we consider 
neutrino scattering). 



dimension |30 . 





(J ~ 



(1) 
(2) 




) 
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while the anti-symmetric piece defines the polarized structure functions gi and g2 

and are the nucleon momentum and spin four-vectors, u = q ■ p, s'^ = —mj^ is our 
choice of normalization, x = Q^/2z/, = —q^ > and m^r is the nucleon mass. 

At the leading twist in the operator product expansion of the two electromagnetic currents 
in Eq. |2l the moments of the structure functions can be factorized, at scale /i, into hard 
perturbative contributions (the Wilson coefficients) and low energy matrix elements of local 
gauge invariant operators. Adopting the notation of js^ . 

2 f dxx-'^F,{x,Q') = 5^cK(/xVQ^^?(/i))^;(^)(/x), 

q=u,d 

[\xx--'F,ix,Q') = J2^t(^''/Q'^9i^^))vi'\^^), 

•^0 q=u,d 

2 f dxx-gr{x,Q^) = ^ eS^i(/xVg', a(f)(/i), 

q=u,d 

2 f dxx-g,M') = $^[ea(^VQ',^7(/^))4^H/^)- 

q=u,d 

- e[%VQ',g{i^))a^^\f,)] (4) 



where cfl, e^fl are the Wilson coefficients and Vn\fi), an^^^\ dn\fi) are the non-perturbative 



matrix elements. At the leading twist Vn {jj) and aW are related to the parton model 
distribution functions {x^)q and {x'"')Aq- 

{x--'), = v'^^ = ^al") (5) 

To determine Vn\lj), cin \ and dn\j^) we need to compute non-perturbatively the follow- 
ing matrix elements: 

I s\Ol^^^,„^Jp, s) = 2v^^\ii) X [p^^p^^ ...p^^ + ...-tr] 

s 

-(P'*l^{'/.iM2...^„}b'*) = ;^^n\f^) X [s<xPmiPm2 ■ ■ - Pm^ + tr] 

{P,s\0fj^^^^^^,„^jp, s) = :^^dl{fi) X [{s^p^, - s^,p^)p^, ■ ■ ■ Pm„ + tr] 

(6) 
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{} implies symmetrization and [] anti-symmetrization of indices. The nucleon states \p,s) 
are normalized so that {p, s\p', s') = {2TiY2E{p)5{p — p')6s,s'- The operators O are 



^w/.2-/.„ = ( 2 ) ^^'^i ■ ■ • Di,„ q - trace, 



= o ^^-^5 D,,--- D,^ q - trace, (7) 



where D — D and D, D are covariant derivatives acting on the right and the left 
respectively. 

In Drell-Yan processes the transversity distribution (x)sn can be measured (for details 
l84l . I85| ) . The relevant matrix element is 

2 

(P' ^l^p.VlM2-Mn}IP' ^) = — ^ [(^PP- - S,Pp)Pf,,P^, ■ ■ ■ PMn + ^'^l 

where 



n 



CpliM2- Mn = ( ^ ) ?750^p^ I^Pi ■ ■ ■ 9 - trace. (8) 



III. LATTICE MATRIX ELEMENTS 

The non-perturbative calculation of the matrix elements described in the previous section 
(structure function moments) proceeds, as do all Euclidean lattice calculations, through the 
computation of nucleon three- and two- point correlation functions, 

C5f(p,t,r) = Y,^-'^{JMt)O{T)UP,0)), (9) 

a,f3 

C2,tm = Ef^^l {Mp,t)Up,o)), (10) 

a,l3 V ^ /a/3 

where J{p, 0) and J{p, t) are interpolating fields with the quantum numbers of the nucleon 
and definite momentum. F"'^ is a Dirac matrix projection operator which is taken as 

r = (11) 



for unpolarized matrix elements, and 

1+74 



^757fc, (^7^4) (12) 
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for polarized matrix elements. 



1+74 



projects out the positive parity part of the baryon 



propagator. For the proton a typical choice for the interpolating field is (36 1 
JM t)=Yl e-'P-^e"'' [uliy^, t)C754(z/2, t)] u,,^{ys, - x)0(i/2 - x)4>{y3 - x) (13) 



x,a,b,c 



with charge conjugation matrix C = 7472, a a spinor index, and a, b, c color indices. The 
functions 0(a;) are smearing functions that are designed to maximize the overlap of the 
interpolating field with the ground state of the nucleon. For the source we used 0(x) = 1 
when X is within a box of size i? ~ 1 fm and zero outside this box. For the sink we took 
a point sink, 0(x) = 6{x). We have optimized the size of the box to maximize the overlap 
of the source to the ground state. This setup works well for zero spatial momentum of the 
proton, and since we studied only this case so far, this was all we needed to do. For non-zero 
momentum this smearing is not optimal; one must resort to other smearing methods such 



3- 



as gauge invariant Gaussian or Wupertal smearing 

In the limit when the Euclidean time separation between all operators is large, t ^ r ^ 0, 
the desired matrix element between ground states dominates the correlation function, 

En{p) + rriN 



C2pt{P,t) 



Zn 



.g -Ei,{p)t 



2E{p) 

Sr^ Uaip,s){p,s\O\p,s')u0{p,s') 



(14) 



where u{p,s) is the nucleon spinor satisfying the Dirac equation, and {0\Ja{p,t)\p, s) = 
\/ZNUa{p, s). Using Eq.HUand Eq.|ni(or Eq.lH)) the desired matrix elements can be extracted 
from the ratio of three point functions to two point functions. In practice we would like to 
achieve the asymptotic behavior of Eq. El with as small as possible t. For that reason the 
smeared interpolating operator J is essential. For more details on the technical aspects of 
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15 



37|, 



the lattice calculation the reader may refer to 

The momentum fraction {x)q carried by each valence quark in the nucleon is computed 
in Euclidean space by inserting into the correlation function the operator 



44 = 9 



74 D4 



^^7fc Dk 



(15) 



D is the lattice covariant derivative, or difference operator. 



with 



= ^ [Uf^i^hi^ + A) - Ulix - fi)q{x - jl)] 



(17) 



C44 belongs to the 3^ representation ^ of the hypercubic group H(4) and does not mix with 
any lower dimensional operators under renormalization The ratio, 



r cT' 







2pt 



mN{x) 



(18) 



then yields the momentum fraction. 

Similarly, the helicity distribution (x)Aq for each valence quark is computed from the 
operator 



O 



5g 



{34} = ^^75 



73 Da +74 -D3 



(19) 



belonging to the 63 representation of H(4). It also does not mix with lower dimensional 
operators under renormalization js^, [s^ . The ratio yields 



R 



{^)Aq 



r (34} 



C, 



2pt 



mN{x) 



Aq- 



(20) 



The lowest moment of the transversity, {l)Aq, related to the tensor charge of the nucleon, 
is computed in Euclidean space using the operator 



34 = g75C^34g, 



(21) 



with cr^j, = |[7^,7i/]. This operator belongs to the 6^ representation of H(4) and does not 



mix with any lower dimensional operators under renormalization 



HQ. 



Again 



R 



Co, 



(1) 



2pt 



Sq- 



(22) 



Finally the twist-3 matrix element di related to gi and g2 is given by the operator 



^5q 



[34] = ^975 



73 Da -74 D3 



(23) 



^ The representations of H(4) are denoted as d^, where d is the dimension of the representation, C is 
the charge conjugation and the subscript distinguishes between different representations of the same 
dimensionahty and charge conjugation. 
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belonging to the 6^ representation of H(4). O^^^-^ is allowed to mix with lower dimensional 
operator O^l if the lattice fermions do not respect chiral symmetry. The mixing coefficient in 
this case is linearly divergent with the inverse lattice spacing and hence a non-perturbative 



14| . The use of domain wall fermions eliminates this problem as first 



subtraction is required 
shown in and in Section IVII di is given by the ratio 



5q 



Rd, = = d,. (24) 

IV. RENORMALIZATION 

The local operators discussed in the previous section arise from an operator product 
expansion and therefore must be renormalized. The renormalized operators defined at scale 
are obtained from lattice-regularized operators defined with lattice spacing a. 

Oiifx) = Zi{fx,a)Oiia) + % (/i, a) O,- (a) , (25) 

where Oj are a set of operators allowed by symmetries to mix among themselves, and dj 
is the dimension of each operator. If mixing with lower dimensional operators occurs, the 
mixing coefficients are power divergent, and hence must be computed non-perturbatively 
to accurately subtract them. The mixing of lattice operators is more complicated than 
that of the continuum operators since not all of the continuum symmetries are respected on 
the lattice. In particular, 0(4) rotational symmetry in Euclidean space is broken down to 
the hypercubic group -ff(4). As a result, an irreducible representation of 0(4) is reducible 
under H{4) and hence mixing of operators that would not occur in the continuum can 
occur on the lattice. For a detailed analysis of the H{4) group representations see 
and references therein. The lattice operators are selected carefully so mixing with lower 
dimensional operators does not occur and thus no power divergences are encountered. In 
general, the breaking of rotational symmetry makes the calculation of higher spin operators 
(and thus higher moments of structure functions) difficult. 

The breaking of chiral symmetry {e.g., by the lattice) results in mixings with lower 
dimensional operators for the dn matrix elements. The problem is avoided by using chiral 
lattice fermions such as domain wall, overlap, or fixed point fermions. 

In order to renormalize the quark bilinear operators studied here, we employ the non- 
perturbative renormalization (NPR) method introduced in This method has been 

8 



shown to work very well in the case of domain wall fermions |41|, |42|. For quark bilinears 
without derivatives we only need to compute a single quark propagator from a point source; 
the Fourier transform then yields the quark propagator S{pa; 0) in momentum space which 
is sufficient to calculated all of the needed vertices. 

S(pa;0) = 5^e-^^-^S(x;0), (26) 

X 

where the lattice momentum is 

2n 

P,. = -r^,. (n^ = 0,±l,±2,...), (27) 

where is the linear lattice size in direction /i. Because the NPR method relies on matrix 
elements of the operators between off-shell quark and gluon states, the calculation proceeds 
in a fixed gauge which, for convenience, is chosen to be the Landau gauge. 

In the case of operators containing derivatives, a single quark propagator is not sufficient 
to build the quark bilinears we need. One way to proceed is the method in where prop- 
agators with momentum sources instead of point sources are prepared in Landau gauge and 
derivative (difference) operators are constructed at the sink point where values at all neigh- 
boring points are available. This method works well when operators with many derivatives 
are needed but is expensive if many values of momentum are desired. 

In our case we are only interested in operators with one derivative which can be con- 
structed from a point-source propagator and a point-split-source propagator with appropri- 
ate gauge links attached at the source. 

S,{pa; 0) = J2 e-'^-^'l [S(x; -A)^/^!-^) - S(x; • (28) 

X 

Following this strategy, we compute the matrix elements of derivative operators with many 
values of momentum by computing four additional quark propagators on each gauge config- 
uration. 
Following 



41|, we define the amputated, bare, vertex function for each operator O, 



and a corresponding projector that enforces the tree-level renormalization condition 

TrVV^ oc 1. (30) 



We choose the following, 

C|44} ^ Vl^'^ = 74P4 -\Y1 

i=l,3 

^{34} ' ^34 ^ = ^75 ^[73^4 + 74^3], (32) 

— n'"' = 75Cr34. (33) 

Defining the renormalized operator at scale /i^ = as OrenifJ') = Zo{fJ',a)0{a), the renor- 
malized vertex takes the form 

V'renip') = {p% (34) 



and after projection, 



Ao.e„(p^) = ^ Aoip') = T^r^^Tr {Vo{v)V''{j?)\ (35) 

From this we extract the required renormalization constant Zq that ensures the same renor- 
malization condition at scale /i as in the free case. 

Ao.en(/i') = ^ Ao(/^') = l- (36) 

Zg 

A. Anomalous dimensions and matching 

In order to compute the renormalization group invariant (RGI) constants [4^, we divide 
out the running of the operator at hand, calculated in continuum perturbation theory. To be 
sensible, this is done at large enough momentum so the operator runs perturbatively. The 
momentum cannot be taken too large, or lattice artifacts will spoil the continuum running 
of the operator. Thus the momentum scale where the RGI constant is defined must satisfy 

Aqcd < < ^- (37) 
In addition, it is convenient to match the RI/MOM regularization used in NPR with the 



MS scheme since the latter is conventionally used in the Wilson coefficient calculation. For 



the derivative operators the matching from the RI/MOM scheme in Landau gauge to MS 
can be done using the conversion factor [l^ 



Gn + Sn-\ 



2(n- 1) iEmPm^m(p) 

10 



Oial) , (38) 



where Cp = 2N quadratic Casimir for the SU{Nc) gauge group, the strong 



couphng constant, 

2 



2 " 1 

- (25„+i - 3 - + - 4 5^ - (25, - - 1 , (39) 



n{n + 1) n + 1 j 



and 



hf,{p)= (^^^^^2■■■^inPf^2■■■P^^n■ (41) 

The conversion factor ^mom depends on the direction of the momentum p and on the coeffi- 
cients c^i,...,^„ because the renormahzation condition breaks 0(4) invariance. The coefficients 
C/xi,...,/i„ are defined using the conventions in In the present case for the momentum 

fraction and hehcity operators we have 

3 

3 



C^,u = ^iii^uA ~ ^ ^tik^uk , (42) 
k=l 

and 

CfMu = 5^35^4. + S^^Sus, (43) 

respectively. 

For matching the tensor bihnear we use results from 



(4320C(3) - 4815) - 1252Tj.nj + (5987 - 3024C(3)) Ca 



(44) 



where ({k) is the Riemann zeta function, Ca = Nc and = 1/2 for the SU{Nc) group. 



In our analysis, we first match to MS and then use the MS running of the operators 
to extract the RGI renormahzation constant. 

Following the conventions of [4^, the continuum two loop running of the operators is 
parametrized by 
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where the anomalous dimension of the operator is 

7o = ^lo ( 



47r/ 



lOi 



(0 
To 

A 
/So' 



(46) 

(47) 
(48) 



and /?o,i are the first two coefficients in the weak couphng expansion of the beta function 



/3(as 
47r 



-/5o 





2 




-47r. 


-A 


.47r. 



(49) 



The two loop running of is given by 

1 



45 



A In In [liV^lcD) 



(50) 



47r /5oln(/iVA|^^) /?o' In^ (/iVA^cfl) 
The values of the (3iS and the 7j's used in this analysis are given in Tables HI and ITTl respec 
tively. For Aqcd we take the quenched value Aqcd = 238 ± 19MeV. 



V. NUMERICAL DETAILS 

We work in the quenched approximation and use domain wall fermions to compute the 
matrix elements described in the previous sections. We use the DBW2 gauge action 0, li^ 
at lattice spacing 1.3 GeV (/3 = 0.87), with lattice size 16^ x 32 and fifth dimension 

Ls = 16. This action has been shown to significantly reduce the explicit chiral symmetry 
breaking of domain wall fermions with finite fifth dimension relative to the Wilson gauge 



action 



This relatively coarse lattice spacing was chosen to give a large physical volume 
(2.4 fm spatial size) and enable calculations with light quark masses to study the chiral 
behavior of the matrix elements. An earlier calculation of the nucleon axial charge showed 



that such a large volume is necessary to avoid significant finite volume errors j2ll| which, based 
on that study, we may expect to be a few percent for pion masses in the range 390MeV to 
850MeV. The residual quark mass for Lg = 16 at this lattice spacing and for the DBW2 
action is mres ~ 0.7 MeV^|, truly negligible compared to the input quark masses in our 
simulation, 0.02 < mr < 0.1, which span a range from about one-quarter to two times the 
strange quark mass|3ll|. 
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To calculate three point correlation functions we use box sources with size ~ 1.2 fm and 



sequential-source propagators with point sinks. For details see j2l|. The source time is 
t = 10 and the source-sink time separation is 10 lattice units, approximately 1.5 fm, which 
provides a sufficiently large time separation to observe clear plateaus. Figures ^0 show 
typical plateaus for the matrix elements studied in this work. Our calculation is done using 
416 independent gauge configurations produced using the overelaxed heatbath algorithm 
described in All statistical errors are jackknife estimates. 

For the calculation of renormalization constants we use 120 lattices and fix to Landau 
gauge using the technique described in Two loop continuum running is used to extract 
the RGI renormalization constants in all cases. In order to eliminate remaining scaling 
violations we fit the data linearly in (ap)^ and define the renormalization constant from the 
intercept (as done in 4l|). The fitting range used is (pa)"^ G [1.2, 1.9]. 



VI. RESULTS AND DISCUSSION 



The nucleon mass has been determined on a subset of this ensemble of lattices 
previously 3l| . For the present analysis, we fit the two-point correlation function to a single 



exponential from time t = 6 to 15. The fitted nucleon mass for each quark mass is given in 
Table llllj they are consistent with those reported in . 

The central question we sought to answer with this study is why lattice results for the 
first moments of the polarized and unpolarized structure functions disagree with fits to 



experimental measurements 



13 



15 



25 



50 



5l|. Preliminary results from this study and 



i3,yQ- 



one using two fiavors of dynamical domain wall fermions have been reported in 
The discrepancy is large 50%), and holds for dynamical as well as quenched calculations. 
A plausible explanation is that the quarks simulated in these past studies, valence and sea, 
have been too heavy. Here we use as light a quark mass that has been simulated to date for 
nucleon structure calculations, roughly one-quarter of the strange quark mass, in order to 
investigate the chiral regime. If the problem is related to the sea quark mass, a resolution 
will have to wait for future dynamical fermion calculations using lighter sea quark masses 
than have been used already^. 



^ The RBC and UKQCD collaborations are embarking on a large scale project this year to generate an 
extensive ensemble of 2+1 flavor domain wall fermion lattices with sea quark masses as light as 1/5 the 
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In Figure El and Table IIVI we present our results for the momentum fraction, of 
the up and down valence quarks. These results do not contain disconnected diagrams. 
Figure IHl shows the isovector matrix element {x)u~d where disconnected diagrams do not 
contribute (for degenerate u, d quarks). The central values are obtained from a constant fit 
over the range 13 < t < 16, based on the plateau in Figure Q (the same range is used for 
all matrix elements in this work). The quark mass dependence is mild and appears linear. 
The renormalization constant corresponding to the momentum fraction is Z'^^^{2 GeV) = 
1.02(10), essentially one, (see Figure IT^ and Table IVlTl) . The renormalized results are all 

nnn 

significantly higher than the value 0.154(3) extracted from the experimental data of 

by Dolgov et.al. Q|, and there is no apparent curvature as that leads us to believe 

that in the chiral limit the two results would agree. Thus, we do not attempt to extrapolate 



to = 0. This is the same behavior witnessed in previous studies; in particular, our MS 
values for {x)q are quite consistent with the recently reported quenched improved Wilson 
fermion, continuum and chiral limit, value in [2^, suggesting good scaling of domain wall 
fermions. 

Figure m shows the first moment of the helicity distribution (x)Ag for the up and down 
quarks (tabulated in Table |Vj). As before, these results do not contain disconnected dia- 
grams, and in Figure |H1 we display the isovector matrix element {x)/\u~Ad where disconnected 
diagrams do not contribute. Again, the quark mass dependence is mild and appears linear, 
the renormalization constant is essentially unity, Z^^^{2 GeV) = 1.02(9) (see Figure HH and 



Table IV11|) , and the renormalized moments 



ie above the value 0.196(4) extracted from the 



Since chiral symmetry requires that the renormalization constants of the momentum 
fraction and the helicity distribution be the same, we can consider the ratio of the bare 
matrix elements in which the renormalization constants and matching factors cancel. A 
.,„,.a. .atio wcUed we.. t.. case of the axia. e...ge. 0. and we saw already that 
the explicit calculations of these constants gave the same result well within statistical errors. 



mm. 



Figu re p shows the ratio together with the value extracted from experiment! 
R B lisjl. As it is argued in [3], it is difficult to estimate the systematic errors associated 
with the experimental extraction of the values of both {x)u-d and {x)Au-Ad but it is almost 



strange quark mass. Among many other quantities, nucleon structure will be studied 
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certain that these systematic errors are smaller than the statistical errors of the lattice 
calculation. For that reason the comparison of lattice results to experiment is not out of 
order. Interestingly, our results are in very good agreement with experiment; a discrepancy is 
not evident at all. Note that the jackknife determination of the ratio is determined relatively 
precisely since the error on the renormalization constants does not enter and the numerator 
and denominator are highly correlated. 

In the axial and vector matrix elements displayed very different finite volume behavior 
as the quark mass was reduced; the axial matrix element, being sensitive to low energy 
physics, decreased drastically in the limit ruf in the small volume study. If similar 
behavior is operative here, the mild quark mass dependence of the matrix elements suggests 
that finite volume effects are small. 

In the quark mass dependence near the chiral limit of {x)u—d but not of (a^)A«— Ad; 
was computed in quenched chiral perturbation theory, so we can not say what the chiral 
perturbation theory prediction is for the ratio. For {x)u~d alone, the chiral perturbation 
theory formula allows a wide range of values at the physical point, depending on the values 
of various unknown low energy constants. The authors of investigate the quark mass 
dependence for several values of the unknown parameters. The most natural choice demon- 
strates a smaller dependence on the quark mass than that predicted by full QCD chiral 
perturbation theory. 

Our results suggest that whatever systematic error causes the discrepancy in the indi- 
vidual moments from experimental expectations appears to mostly cancel in their ratio. 
Although our calculation is quenched, it is instructive to look at the full QCD chiral per- 
turbation theory formulas found in j5l|, |52| • 



{x)u-d = C 

{x)Au-Ad = C 



(51) 
(52) 



where the normalization is such that the physical pion decay constant is /tt = 93 MeV, 
C and C are unknown constants, and e(/x^) and e(/i^) are counter terms evaluated at the 
renormalization scale /x. In Figure El we plot the above formulas for /i = 1 GeV. The 
unknown constants C and C are chosen so that the formulae reproduce the experimental 
result at the physical pion mass point, and the counter terms are set to zero in order to 
isolate the effect of the chiral logarithm. We see that there is a strong dependence on the 
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pion mass. As a result, when m,r ~ 400 MeV, the momentum fraction is roughly 50% larger 
than at the physical point, while the first moment of the helicity is about 30% larger. This 
indicates that the discrepancy between the lattice data and experiment may be due to the 
unphysically large masses used in current lattice simulations. On the other hand, the large 
size of the one-loop perturbative corrections suggests that chiral perturbation theory at this 
order is unreliable. The effects of the counter terms and higher order contributions may also 
be large and could tend to cancel the large one-loop corrections. Nevertheless, it is worth 
noting that our quenched results at pion masses of about 400 MeV differ from experiment by 
roughly the same amount as one-loop chiral perturbation theory in full QCD predicts. The 
observation of the large one-loop corrections was first made by the authors of [26, 53, 54, 5^. 

On the other hand the ratio of the momentum fraction to the first moment of the helicity 
(Figure IT7|) is a milder function of the quark mass. The difference of the experimental result 
from the value at = 400 MeV is about 10%, not wildly inconsistent with the lattice 
result shown in Figure 01 

The conclusion from this discussion is that the discrepancy between lattice and experi- 
mental results for the quark momentum fraction and helicity distribution is most likely due 
to strong mass dependence in these functions and will be resolved by pushing lattice simula- 
tions further into the light quark mass region. Since dynamical results so far are similar to 
quenched 3; , this probably holds in that case too. It may also be of interest to obtain 
the next higher order results in chiral perturbation theory for these quantities, though this 
seems a more daunting task. 

We now turn to another interesting feature in our calculation, the twist-3 matrix element 
di- Although it is not measurable in deep inelastic scattering of electrons on protons, it 
serves as an example of what can be expected for the (i„ matrix elements. As discussed 
in Section HVl the operator used to calculate di, C^j'^j,], mixes with the lower dimensional 
operator O'^l when chiral symmetry is explicitly broken. With domain wall fermions, unlike 
Wilson fermions, chiral symmetry is not broken {rrires is small enough to ignore in this 
study), so power divergent mixing should not occur. Results are summarized in Table IVll 
Figure ITUl shows, unlike the Wilson fermion result this matrix element is small in the 
chiral limit. In fact, the power divergence in the Wilson fermion case switches the sign of 
the u and d quark contributions. Using Wilson fermions, the QCDSF collaboration found 
similar results for ^2 after a non-perturbative subtraction of the power divergence j38t|. These 
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results confirm our expectations that the domain wall fermion formulation avoids the power 
divergence present for Wilson fermions. Because the value of di computed here and QCDSF's 
value of d2 appear small in the chiral limit, we conclude that the Wandzura-Wilczek relation 
between moments of gi and g2 |52t], which assumes vanishing dn, is at least approximately 
true. This relation is not obvious in a confining theory js?! . 

Finally, we have computed the first moment of the transversity distribution, (l)^^, for 
up and down quarks (Figure HTJ and the isovector combination (l)^^-^^ f Figure IT^. The 
transversity is an important target of the RHIC spin program at Brookhaven National Labo- 
ratory (see and references therein). In both cases the quark mass dependence is mild 
and appears to be linear. The renormalization constant (Figure^! and Table EH} is Z^^{2 
GeV) = 0.872(11). Naive linear extrapolation to the chiral limit yields (l)^^-^^ = 1.193(30) 



for /i = 2 GeV in the MS scheme. This result contains an unknown systematic error from 
the chiral extrapolation which may be small given the mass dependence exhibited in Fig- 
ure El We note that the calculation of the flavor non-singlet tensor charge is similar to q a 
which also exhibited mild mass dependence and whose value agrees well with experiment j2l| . 
We caution the reader that the momentum fraction and helicity distribution also exhibited 
mild mass dependence but are known to disagree with fits to the experimental data. Our 
value is consistent with another recent calculation jo^. 

VII. CONCLUSIONS 

We have reported on a quenched calculation of the first moments of the polarized and 
unpolarized structure functions of the nucleon. We have used domain wall fermions in a 
relatively large physical volume, ~ (2.4 fm)^. The large volume is important to minimize 
finite volume errors in nucleon matrix element calculations as shown earlier in the calculation 
of the axial charge, qa I^^- The chiral symmetry of domain wall fermions may also be 
useful in unraveling the mystery of the large discrepancy between lattice calculations and 



experiment for the moments of structure functions. While our results 



"or the individual 



moments show no evidence of the chiral logarithm predicted in |50l l5lL l52l 61 1 for pion 
masses as low as 390 MeV, the ratio of the momentum fraction to the first moment of the 
helicity distribution is in very good agreement with the experimental value (Figure Ej). We 
note that the ratio is computed on the lattice more accurately than the individual moments. 
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This agreement, taken together with chiral perturbation theory calculations, leads us to 
conclude that the discrepancy between lattice calculations and fits to experiment is most 
likely due to strong mass dependence of these functions in the light quark regime. Thus, 
ultimately the difference will be resolved as lattice calculations push into this regime. 

A recent large scale, detailed, investi gati on of the discrepancy in low moments of nu- 
cleon structure functions was reported in 25| ; using improved Wilson fermions, the authors 
performed continuum and chiral extrapolations but did not resolve the problem. The agree- 
ment with our results suggests that scaling violations for domain wall fermions are mild (the 
coarsest lattice spacing used in [2^ was about 0.1 fm compared to 0.15 fm used here). 

Our calculation of the di matrix element (Figure irUj) indicates that the domain wall 
formalism eliminates power divergent mixing in this class of matrix elements. This suggests 
that a precision calculation of d2 with domain wall (or other chiral) fermions is possible 
and should be undertaken. Since the mixing with lower dimensional operators induced by 
explicit chiral symmetry breaking is linearly divergent with a~^, care must be taken to 
minimize rrires- 

Finally, we have computed the first moment of the transversity distribution, or tensor 
charge, which will be measured at Brookhaven National Laboratory as part of the RHIC 



spin program. We find {l)5u-5d = 1.193(30) at /i = 2 GeV in the MS scheme, with unknown 
systematic error stemming from the linear chiral extrapolation. The mild quark mass depen- 
dence in the tensor charge suggests this systematic error is small, as does the similarity to 
the calculation of qa- On the other hand, such mild mass dependence is also observed in the 
momentum fraction and helicity distribution, hence further study of the chiral extrapolation 
of this observable is needed in order to obtain a reliable estimate of the systematic error. 

An unknown systematic error due to quenching exists in our results. Two flavor calcu- 
lations using Wilson fermions with relatively heavy quark masses exist jisl Ig^ . we have 
begun two flavor domain wall fermion calculations Q|, and 2+1 flavor domain 

wall fermion calculations are just beginning. These studies begin to address the quenching 
error. In addition, systematic uncertainties due to continuum, chiral and infinite volume 
extrapolations must be addressed as in the extensive quenched study using improved Wil- 
son fermions reported in 25^. This study suggests domain wall fermions should facilitate 
these extrapolations. Recent results from chiral perturbation theory j64l. 16 5| may also prove 
useful. 
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TABLE I: The one- and (avo-1oo|:) oxpausion c(K\[ricicii(.s of tho fiuouchc-d QCD :3 function. 
Pi quenched (n/ = 0) value 

Po 11 

Pi 102 



TABLE IL Anomalous dimensions of quark bilinear operators used in this paper in quenched 
(n/ = 0) QCD. 

Operator 7o 7i 

OI4 64/9 96.69 

Oj^4j 64/9 96.69 

O3I 8/3 724/9 



TABLE III: The nucleon mass in lattice units. Values are from fully covariant, single exponential, 
fits to the two-point correlation functions in the range 6 < t < 15. Errors are statistical only. 



ruf mjv (error) (dof) 

0.02 0.856 (10) 5.9 (8) 

0.04 0.967 (6) 1.9 (8) 

0.06 1.064 (5) 3.5 (8) 

0.08 1.155 (4) 5.1 (8) 

0.10 1.241 (4) 6.4 (8) 
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TABLE IV: Lowest moment of unpolarized structure functions for u and d quarks and the flavor 
non-singlet combination. Errors are statistical only. Values are not renormalized. 



rrif {x)u {x)d {x)u-d 

0.020 0.446(18) 0.184(11) 0.262(16) 

0.040 0.446(7) 0.191(4) 0.255(6) 

0.060 0.457(5) 0.2012(28) 0.256(3) 

0.080 0.464(4) 0.2092(21) 0.2548(25) 

0.100 0.471(3) 0.2160(18) 0.2549(21) 



TABLE V: Lowest moment of the polarized structure fucntions. Errors are statistical only, u and 
d valence quark contributions. Au — Ad denotes the flavor non-singlet combination. Errors are 

statistical only. Values arc not rcnonualizccl. 



•nif {x)au {x)Ad {x)Au-Ad 

0.020 0.271(21) -0.060(13) 0.331(23) 

0.040 0.261(7) -0.069(4) 0.330(8) 

0.060 0.260(4) -0.0671(23) 0.328(5) 

0.080 0.263(3) -0.0664(16) 0.330(4) 

0.100 0.2667(26) -0.0667(12) 0.3335(29) 



TABLE VL The twist-3 matrix element di. u and d valence quark contributions. Errors are 
statistical only. Values are not renormalized. 

rrif diu did 

0.020 0.045(21) -0.029(14) 

0.040 0.075(8) -0.031(5) 

0.060 0.109(5) -0.0354(29) 

0.080 0.140(4) -0.0405(20) 

0.100 0.167(3) -0.0457(15) 
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TABLE VII: The MS renormalization constants at = 2 GeV. Errors are statistical only. 

^(.), 1-02(9) 
%>A, 1-02(10) 
Zt 0.872(11) 
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FIG. 3: The three point correlation function for the operator C[34] and bare quark mass 0.040. 
Octagons are the up quark contribution and diamonds are the down quark contribution. 
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FIG. 4: The three point correlation function for the operator O^l and bare quark mass 0.04. 
Octagons are the up quark contribution and diamonds are the down quark contribution. 
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FIG. 5: The bare momentum fraction, m.,^ is the pseudoscalar mass. Octagons are the up quark 
contributions and diamonds are the down quark contributions. Disconnected diagrams are not 
included. 
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FIG. 6: The bare flavor non-singlet momentum fraction. 
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FIG. 7: The bare helicity distribution. Octagons are the up quark contributions and diamonds are 
the down quark contributions. Disconnected diagrams are not included. 
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FIG. 8: The bare flavor non-singlet helicity distribution. 
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FIG. 9: The ratio of the flavor non-singlet momentum fraction to the helicity distribution (oc- 
tagons). The experimental expectation is marked by the burst symbol. 
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FIG. 10: The bare di. Octagons are the up quark contributions and diamonds are the down quark 
contributions. Disconnected diagrams are not included. 
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FIG. 11: The bare transversity. Octagons are the up quark contributions and diamonds are the 
down quark contributions. Disconnected diagrams are not included. 
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FIG. 12: The bare flavor non-singlet transversity. 
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FIG. 13: Renormalization constant for the momentum fraction. The diamonds are the renormal- 
ization group invariant points. The octagons are the raw data. 
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FIG. 14: Renormalization constant for the helicity distribution. The diamonds are the renormal- 
ization group invariant points. The octagons are the raw data. 
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FIG. 15: Renormalization constant for the transversity. The diamonds are the renormahzation 
group invariant points. The octagons are the raw data. 
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FIG. 16: The leading chiral logarithm dependence for the first moment of the hehcity and the 
momentum fraction. The curves are normalized so that at the physical point the experimental 
result is recovered. 
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FIG. 17: The leading chiral logarithm dependence for the ratio of the momentum fraction to the 
first moment of the helicity. The curve is normalized so that at the physical point the experimental 
result is recovered. 
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